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We consider the recent proposal [l| for the extraction of the graviton propagator from the spinfoam 
formalism. We propose a new ansatz for the boundary state, using which we can write the prop- 
agator as an integral over SU(2). The perturbative expansion in the Planck length can be recast 
into the saddle point expansion of this integral. We compute the leading order and recover the 
behavior expected from low-energy physics. In particular, we prove that the degenerate spinfoam 
configurations are suppressed. 

I. INTRODUCTION 

Recently [H, [1; H, S] j a proposal has appeared for the computation of the 2-point function of quantum gravity within 
the spinfoam formalism [6| , a candidate covariant approach to a non-perturbative quantisation of General Relativity 
(GR). This proposal offers a possibility to study the semiclassical limit of spinfoams and define the perturbative 
expansion in the Planck length £p, arguably the major open question within this approach. This is very interesting, 
for two different reasons: on the one hand, to check the correctness of the low-energy limit of spinfoams; on the other 
hand, to address the possibility of curing the non-renormalisability of the conventional perturbative expansion via 
background-independent methods. 

In in (see also [sj), the proposal was applied to the Barrett-Crane (BC) spinfoam model for four dimensional (4d) 
Riemannian GR, and it was shown that the 2-point function, or graviton propagator, indeed had the correct l/\x — y\'^ 
leading order behavior in the large scale limit. However, there are some assumptions behind this result that deserve 
a more careful treatment. In particular, the BC model is plagued by degenerate configurations. This has so far cast 
some doubts on the viability of the proposal, especially because there is still no support from the numerical analysis, 
due to the high complexity of the vertex amplitude. 

In this work, we prove that the degenerate configurations do not affect the leading order. In doing so, we modify the 
ansatz for the boundary state, following the 3d investigation appeared in For the 4-simplex spinfoam contribution, 
the new boundary state allows us to write the 2-point function as an integral over SU(2). This is our first result. In 
the large j limit, we evaluate this integral in a saddle point approximation, proving the l/\x — behavior of the 
leading order. The degenerate configurations can be correctly neglected because they correspond to saddle points 
which are not absolute minima. This is our second result. Thanks to the integral expression, which involves no sums, 
the numerical analysis is strongly simplified. Indeed, the leading order result is fully supported by numerical analysis 

Concerning the full perturbative expansion in £p, this can be computed from higher orders of the saddle point 
approximation as well as from the contribution of other spinfoams. We do not attempt here their evaluation, which 
is rather complicate. Notice however that higher orders have been studied in a 3d toy model 0]; remarkably, 
interesting modifications to the conventional expansion arise, due to the microscopical quantum geometry described 
by spinfoams. This makes the whole approach very interesting, and pushes towards the calculation of these corrections 
in the 4d case. We expect that having recast the £p expansion into the saddle point expansion will help future work. 

This paper is organised as follows. In the next Section, we briefly recall properties of the BC vertex amplitude 
which we need in the rest of the paper. In Section IIIIl we review the construction of the 2-point function and the 
results of Q. In Section HVl we introduce the new boundary state and discuss its properties. In Section |Vl we show 
how the new boundary state allows to write the 2-point function as an integral over SU(2). In Section IVIl we study 



* elivine@pcrimetGrinstitute.ca 
t sspezialeOperimeterinstitute.ca 



2 



the saddle point approximation to this integral, and show that the leading order behaves as expected from low-energy 
physics. In the final Section, we discuss possible further developments. 
Throughout the paper, we use units £p = 1. 

II. SPINFOAM AMPLITUDES: EVALUATION OF RELATIVISTIC SPIN NETWORKS 

In this Section we recall basic facts about the spinfoam amplitudes, the quantities encoding the dynamics of 
quantum gravity. This gives us the opportunity to fix our notation and to pinpoint features which will be crucial in 
the following. We consider here the spinfoam amplitude for the BC model of 4d Riemannian quantum gravity; this 
is given by the evaluation of a relativistic spin network 0] with group Spin(4), namely the double cover of SO (4). 
The relativistic spin network is defined by a graph T together with the assignment of a Spin(4) group element d 
to each node n of F and a simple irreducible representation (irrep) J; to each link I oi T. Using the homomorphism 
Spin(4) = SU(2)i x SU(2)fl, the irreps of Spin(4) are labelled by two half-integers, say {j,k), corresponding to the 
irreps of the two SU(2) sectors; then the simple representations are such that they induce the same SU(2) representation 
in the left and right sectors, namely J; = Furthermore, let us note here that the scalar Casimir of Spin(4) 

satisfies C|pi„(4)(j, k) = i [C^{j) + C^{k)], where C^{j) is the SU(2) Casimir. 

Using the above homomorphism, each group element decomposes as the product of two left and right rotations 
G = gLgR, and the evaluation reads Q, 

JSpin(4) „ , 

where s{l) and t{l) denote the source and target node of the link I. The kernel JCj{G) is the matrix element of G on 
the SU(2)-invariant vector \ J, 0) in the J representation. Here SU(2) is the diagonal rotation group, corresponding to 
the subgroup of 3d rotations. We conveniently parametrise SU(2) group elements as 

(7(0, n) — cos (f>l + i sin (jjh ■ a, G [0, tt] . (2) 

Consequently, the characters are given by Xj(ff) = ^^^^ Haar measure is dg = sin^ (j)(Pn{h) d(j). 

The invariant vector is easily expressed in term of left/right components: 



I J, 



0) = ^^(-ir™b-,m)i|.7, -m)fl, = J-Y^\j^m)LR.{3M, (3) 



' di — ' ' \/ d^ 



where dj = 2j + 1 is the dimension of the SU(2) representation of spin j. Then it is straightforward to realize that 
the SU(2) invariant kernel ICj is simply the SU(2) character: 

(J,0|G|J,0) = (J,0|,9l.9h|J,0) = 1^0-,m|5L5j,i|j,m) = ^XjigLg],'). (4) 

Finally, using the invariance of the Haar measure dG ~ dgLdgs, under left and right multiplication, it is easy to prove 
that the relativistic spin network evaluation is actually a 3d object regarding only integrals over SU(2): 



£v = 



/ X{'^^9nX{^X]A9s(l)9m): (5) 

"'SU(2) „ , «7"! 



In particular, the vertex amplitude for the BC model is obtained for F given by a 4-simplex, and this gives the {lOj} 
symbol for the recoupling theory of SU(2), 

{10j}=/ X{d9nX{^Xni9:^i^9m)- (6) 

JSU(2) „ , ji 
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Let us point out that the {lOj} symbol is defined up to a normalisation. This creates an ambiguity in the definition 
of the BC model. The standard normalisation found in the literature is ([6]) without the \/dj factors. Notice also that 
another natural choice for the kernel is cPjlCj = djXj, instead of Xj/dj, as it normalises the convolution product: 

/ dG {d]K.j{HG-^)) {d]K.j{GK-^)) = d]lCj{HK-^). 

JSpin(4) 

These normalisation issues do not modify the computations below, namely the leading order of the graviton cor- 
relation computed using a single 4-simplex. We nevertheless mention them as they will most likely affect the higher 
order corrections and become essential when considering configurations with many 4-simplices. 



Geometrical interpretation of the {lOj} symbol 

The {lOj} symbol admits a geometrical interpretation, associated to the structure of a 4-simplex, which will be 
important in the following. The key fact is that it can be written as an integral over ten SU(2) angles (t>i e [0,7r], 

{m = j d^i[cj>l]\{^X3Ml). (7) 

where the measure takes into account that the angles come from the vertex group elements g„ through the relation, 
cos(j)i = 5tr(p~^J^g((;)). Let us introduce the notation (/J) for I linking the nodes / and J, such that 4>ij = (jji, with the 
convention 0// = 0. The above relation imposes a constraint that can be written as the vanishing of the determinant 
of the 5x5 Gram matrix Gjj = cos^/j, as shown in (l^ : 

MM = lldq^i sin(/)z(5(detG/j). (8) 

This constraint has a clear geometrical interpretation: it says that the angles are the dihedral angles of a certain 
4-simplex. Indeed, notice that the spin network induces a dual triangulation which is also a 4-simplex, with tetrahedra 
dual to the nodes and triangles dual to the links (see FiglT]). Then the constraint can be translated into the Schlafli 
identity d(pi = 0, where Ai is the area of the triangle (dual to the link I) of the geometric 4-simplex, and 

(pi its dihedral angle. The areas can be written as derivatives of the Gram matrix, Ai = k ^^'^g^''' , where k is a 
proportionality constant, related to the 4-volume of the simplex. For more details see the Appendix. 

The original group integration in ^ is over 5 copies of the 3-sphere ^ Spin(4)/SU(2) and the 10 angles are 
easily related to the 5 original 4d unit vectors Nj G S^, via cos^/j = Nj ■ Nj. 

For later use, let us consider the equilateral case when all the dihedral angles are equal, (jji — ip. Then the constraint 
simply reads 

det Gij{(j)i v?) = (1 - cos ip)^ + A cos ip), (9) 

whose roots are given by the completely degenerate 4-simplex ip = 0, and the equilateral 4-simplex ip = 9 — 
arccos(— 1/4). 

To summarise, two important features emerge from the above discussion: 

• The relevant group for 4d Riemannian quantum gravity without matter is SU(2). 

• The angles (f> parametrising SU(2) group elements can be thought of as dihedral angles between the tetrahedra 
dual to the nodes of the boundary spin network entering the spinfoam amplitude. 

Equipped with these considerations, we now proceed to describe how to construct the 2-point function. 
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FIG. 1: The 4-simplex (or pentahedral) boundary spin network. We label the nodes / = 1, . . . 5. In the dual picture, they are 
in correspondence with tetrahedra of the boundary triangulation. Two of them are represented. The links /J, on the other 
hand, are dual to triangles. Consider for instance the link 45: this is dual to the triangle shared by the tetrahedra 4 and 5. 
Associated with the link 45 is the dihedral angle 045 between the tetrahedra 4 and 5. 



In the conventional quantum field theory framework one expands the metric tensor around a given background, 
typically the Minkowski flat spacetime rj^^, writing g^^, = 77^^ + Iph^^] then the 2-point function. 



is evaluated in perturbation theory. The leading order, coming from the quadratic term in the action, goes as l/\x~y\'^, 
namely as in momentum space. The higher orders in the action give self-energy corrections. 

In the spinfoam formalism, the 2-point function (fTO|) can be studied looking at the correlations between fluctuations 
of geometrical quantities. The fluctuations are defined with respect to a boundary metric which is encoded in a spin 
network state s. The points x and y in (|10p can be identified with nodes of s, and the directions ^ can be defined 
using the links of s. In particular, since each link of s can be thought of as dual to a triangle in a triangulation of the 
boundary, it is convenient to project pO[) along the normals to these triangles. Namely, using the boundary metric, 
we contract the indices of (jlOp with four normal vectors, two belonging to the (3d region dual to) the node x and 
two to the node y. For each node, we should distinguish the case when we take twice the same normal, from the case 
when we consider two different normals. In the first case, the contraction with /i^^ gives (the fluctuation of) the area 
of the triangle. In the second case, of the dihedral angle between the two triangles. 

In the quantum theory, both areas and dihedral angles are represented by quantum operators. The area operator has 
spectrum given by twice the square root of the SU(2) Casimir operator,^ A = 2 C{j) = 2j + 1. Analogously, the angle 
operator can be expressed in terms of SU(2) Casimir operators, via the recoupling theorem (see for instance 0, B])- 
This is how the spin variables entering the amplitudes are related to geometrical quantities. Since the fundamental 
variables of the theory are triangle areas and dihedral angles, we compute the 2-point function as the correlation 
between fluctuations of these quantities around a flat background. Choosing a flat background for the boundary 
geometry allows us to introduce a (spin) scale parameter, which we call Jq. We do so by describing the triangulation 
of the flat boundary in terms of equilateral triangles with areas ^0 = 2C(jo). The parameter jo can be used to 
measure the physical distance between two points on the boundary with respect to the flat background metric. 

In the following, we consider only the correlation between the areas, leaving a discussion of the other cases for the 
conclusive Section. Introducing the normal vectors and to the triangles a and b, the area correlations can be 
compared with the projections n'^^n'^n^n'^ Wfii,pa-{x, y) of the continuum 2-point function. Notice that if the 4-simplex 
is equilateral, there are only three independent projections these correspond to the three cases when the triangles 
a and b are the same, when they share a side, or when the share only a point. Following the area correlations are 



^ The factor 2 between the triangle area and the Casimir can be understood as follows. Given the two edge vectors , of the triangle, the 
Spin(4) Casimir can be identified with the norm of the bivector B^'^ = | (a^b^ — b^a"^), namely with the quantity ^ |ap |f)p sin^ 6 = 2A^ . 
Therefore 2A^ = C'sp;n(4) (ii i) = using the homomorphism between Spin(4) and SU(2) introduced above. As for the spectrum, 

let us recall that the Casimir operator is always defined up to an additive constant. Usually one takes this constant to zero, so that its 
spectrum reads C^{j) = j{j + 1). On the other hand, here we take a shift of 1/4, so that C'^{j) = [j + 5)^, to match the results on the 
asymptotics of the {lOj} symbol reported below. 



III. THE 2-POINT FUNCTION 



W^upa{x,y) = {0\T{h^4x)hp4y)}\0) 



(10) 
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given by 



WabUo) = E ^(^-) ^^^b) A[S], (11) 



where J\f = 4*0 [s] Als] is the normaUsation, and the factor 1/jq conies from normahsing the projections along the 
normals. 

Let us briefly explain this formula, referring to Q for a more complete description. The graviton propagator is 
expressed as a sum over all possible boundary spin networks. This includes a sum over all the the graphs and a sum 
over all possible assignments of spins to the links of the graphs. 

The boundary state ^'o[s] should represent a coherent semiclassical state of the boundary geometry In particular, 
we require that the relative uncertainties of the geometry on this state vanish in the large jo limit, namely 

(I'olAjil^o) (^o|A0i|^o) 

(*o|jd*o) " ' (*o|0d*o) ^ ^ ' 

for all links I in the boundary spin network. 

The quantity h.{ja) := C^{ja) — C^ijo) represents the field insertion h^j_^{x) in a Coulomb-like gauge-fixing We 
see that it is (one fourth) the fluctuation of the squared area. 

The amplitude A[s\ is given by the spinfoam model chosen, and it is a sum over all spinfoams a whose boundary 
d<T is given by the spin network s. Working with the Barrett-Crane model, we have 

= E E 114 n ^^(j/) n (^e(.7/))^ X{{m, (13) 

gjLj, is f eecr\d(T e£da v 

where v, f are respectively the vertices and faces of the spinfoam, and jf = ji whenever a face / intersects the 
boundary forming a link I. The vertex amplitude is given by the {lOj} symbol introduced in the previous Section. 
Notice that we left unspecified the edge amplitude Ae{jf). This reflects the normalisation ambiguity of the {lOj} 
symbol, and different choices lead to different versions of the BC model^. However we expect changes in the edge 
amplitudes not to affect the leading order of the graviton correlations Wat, in analogy with the 3d case [|. They will 
nevertheless affect the higher order corrections, allowing to distinguish and discriminate the different choices of A^. 
Below we will choose a particular edge amplitude which simplifies the calculations of the leading order. 

The expression (|lip is the full graviton propagator. This can be evaluated in perturbation theory. The parameter 
of the conventional perturbative expansion is the Planck length £p. By dimensional analysis, it appears (squared) in 
front of jo- We are then led to consider the limit £p i-^ 0, j i— > oo such that £pj is constant, to study the semiclassical 
behavior of the 2-point function. This idea is supported by the asymptotic behavior of the {lOj} symbol for large 
spins, which we recall here. 

In the homogeneous large spin limit, namely when ji = Nki and N i-^ oo, the integral ^ defining the {lOj} symbol 
can be evaluated with a saddle point approximation, leading to the asymptotics [ill , [l^ 

{lOj} ^ E ^(^) (^r[J'] + ^r^) + (14) 

r 

Here S'rJj;] = "^i^jiOi is the Regge action associated with the 4-simplex, dual to the foam vertex, with triangle 
areas Ai = dj,; P{t) and k^- are factors depending on the combinatorial structure of the 4-simplex, and D{ji) is a 
contribution coming from degenerate configurations of the 4-simplex. In principle, the emergence of the Regge action 



^ In the literature, the edge (or tetrahedron) amphtude Ae{ji, j2, ji, ji) is usually the product of some particular powers of dj-^dj^dj^dj 

3ace and is givi 

4 

n Xj^(9)■ 



SU{2) 
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supports the idea that the large j limit can be used to study the semiclassical properties of the theory. However, 
the degenerate term D{ji) dominates strongly the asymptotics, thus hiding the physically interesting first term of 
(fTi)) . Indeed, the numerical analysis in [l^ could not confirm (HH). This fact has raised doubts over the BC model. 
Nevertheless, it was suggested in [l[ that the sick term D{ji) is in fact negligible in the computation of pT|) . because 
it does not match the boundary data induced by [.7;]. Before showing how this happens, let us discuss the structure 
of the perturbative expansion. 

Using the approximation ()14p in we compute the 2-point function within a particular Regge path integral 

formulation of quantum gravity. The measure for this path integral is induced by the spinfoam formalism. On the 
other hand, the Regge action is a discretised version of GR, and it can be expanded around the flat background, 
as in the continuum. The quadratic term in the action is responsible for the behavior of the leading order of the 
2-point function. Higher order terms in the action give corrections in powers of ^p. However, because (|14p is only an 
approximation, we expect to have higher order corrections to the 2-point function which are not described by Regge 
calculus. Indeed, this expectation is confirmed in 3d where an interesting structure of the corrections emerges. 

This is not the end of the story. For any s in the amplitude (jl3[) sums over all spin foams interpolating s. 

Using as in [H the group field theory generated BC model, each foam is weighted by a factor A^, where V is the total 
number of vertices in a, and A is the dimensionless coupling constant of the group field theory. Assuming that the 
perturbative expansion in A is well defined, this suggests that from the group field theory point of view, the dominant 
contributions come from the simplest foams. Indeed, at first order in A we have a single contribution to (|lip (see 0]), 
whose boundary is the 4-simplex spin network. 

Understanding the precise interplay between the (.-p expansion and the A expansion is beyond the scope of this 
paper. In the rest of this paper, we focus on this single contribution. See the conclusions for more comments. 



The Gaussian boundary state and the leading order of the perturbative expansion 

Let us recall here the contribution to the 2-point function coming from a single 4-simplex, as computed in [l[ . When 
we restrict the boundary s to be only the pentahedral spin network, (jlip reduces to 

W^afcOo) = *0[-^''] ^(^-) ^O'^)- (15) 



To explicitly define A[ji\, we need to make a choice for the edge amplitudes Ae. The amplitude A[ji\ will then involve 



the face factors and the square-root of the boundary edge amplitudes (Ag)^/^. We choose Ae(j/) = (D/'^j/) ^ 



Taking into account the five tetrahedra of the boundary of the 4-simplex, we have Ileeao- \/^e(i/) = H/ lie '^j/^^ = 
Y\f dj^ and thus: 

^[^■'1 = n 4 n ^ { 10^} = n {lo^} = / n ^an n iaj'^gm ) . (16) 

To compute the free propagator, namely the leading order in the large jo expansion of ()lip . it is sufficient to consider 
a Gaussian state (see discussion in Q) peaking the (discrete variables representing the) boundary geometry (encoded 
in the spin network labels ji) around a given semiclassical configuration. Choosing the equilateral configuration as the 
background geometry of the boundary 4-simplex, the data encoded by ^'o[j;] are the value 2jo of the area of the ten 
triangles, representing the intrinsic curvature, and the value 9 = arccos(— i) of the ten dihedral angles, representing 
the extrinsic curvature. In particular in the following ansatz was taken: 



M/obi] = cxp I E ^-ji' + + 1)| : (17) 



where 5ji — j/ — jo , and an is a (D-valued (non diagonal) matrix that can be fixed by comparing the leading order of 
(fTT|) with (a suitable discretisation of) the conventional free propagator Q. Using ini), we have '^j^ = -jl a 



(A0i) _ 1 



so that is satisfied in the jo ^ oo limit. 
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It was assumed that the phase of (jl7p suppresses the term D{ji). Then, using 

(HI and (HZl), in Q it was showed 

that the leading order of (jl5[) . in the large spin limit, is 

Wr^—. (18) 
Jo 

This shows that the leading order of the components of the 2-point function behave as expected: rescaling the 
boundary geometry where the two points lie, Wab scales as l/\x — yp. 

However this evaluation relies on a number of assumptions. In particular the use of Ijl4p . and the suppression of the 
degenerate configurations term D{ji). It would thus be important to confirm the analytic result (|18p with numerical 
simulations of (|15p . but this has not been achieved so far [lj|, due to the high complexity of the sum. In the rest of 
this paper we address this issue, and we prove (jlSp in a way that has been verified numerically In particular, 
we do not use (|14p . and we prove that the degenerate configurations plaguing the BC model are indeed suppressed. 
First of all, notice that using the integral representation of the spinfoam amplitude, and swapping the integrals 
with the sums, we can rewrite (|15l) as 

Wab=^^ fH dgn J2 n (9;^l)9m) *0[JZ] Hja) Hjb). (19) 

The key idea is to use a boundary state \E'o[j/]; so that we can perform exactly the sums in (p9)) . To do so, recall 
that the kernel is nothing but the SU(2) character. Then, to be able to perform the sum, it is sufficient to have a 
state V^obi] with a well-defined Fourier transform. This is what we do in the next Section. 



IV. THE NEW BOUNDARY STATE 



As we showed above the relevant group for 4d Riemannian GR without matter is simply SU(2); indeed, pT|) uses 
only SU(2) spins. This has interesting consequences, as SU(2) is also the relevant group for 3d Riemannian GR. 
Therefore we can apply to the Barrett-Crane model in 4d the same techniques developed to study the graviton 
propagator in the Ponzano-Regge model for 3d quantum gravity in In particular, we consider the following new 
ansatz for the boundary state, 

«'o(ji) = n^o(^')' (20) 
I 

Mj) = V"'^"^ '^+-' ^±liiii cos(d,-^). (21) 

V^0(li)-/2,„ + l(^) 

Here the /„(z) are the modified Bessel functions of the first kind, and a E 1R+ is a free parameter.'^ 

Notice that this new boundary state factorises in link contributions and is real, differently from p7|) . As shown in 
0, the Jo I— > oo limit of (|2ip behaves as a Gaussian peaked around jo, 

Mj) - (/^exp{-^(j - jof} cosidjO), (22) 

and thus ((20|) satisfies the semiclassical requirements ([T2|) . Thanks to the asymptotic behavior ((22|) , a boundary state 
constructed from pT|) can be used to reproduce the leading order using the same procedure outlined in the 
previous Section. 



^ With respect to [J, there is an important difference in the phase, which is given by djB here and by dj6/2 in the 3d case. The reason 
for this Ues in the asymptotics I I14I I. which reproduce the Regge action with areas A = dj, whereas the asymptotics for the 3d spinfoam 
amplitude reproduce the Regge action with lengths £ = dj/2. 
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However, the exact form (|2ip opens the way to a new treatment, which aUows us give a check of the procedure of 
The key property of pT|) which we need in the following is its Fourier transform (see [1]), 



2N sin( 



where 



TV : 



a 



(23) 



(24) 



is the normalisation with respect to the Haar measure. From the Gaussian shape of (|23p . it is clear that this new 
boundary state peaks the SU(2) class angle (/) around the values rjO modulus tt. For later purpose, we introduce a sign 
variable cr = ± and write the peaks as 



(j> = r]9 + 



1 



(25) 



As we show below, the different values of the signs 77 and a give the same contribution to the 2-point function, thus 
summing over them simply adds a trivial redundancy. 

The fact that (|2ip admits a simple Fourier transform is at the heart of our construction of the 2-point function as 
an integral over SU(2). In particular, in the following we also need the convolution product between ipo and the field 
insertion h, ipo o h. This can be easily evaluated, using 



1 



■ dft, sin^ (h dck 



V'X,(</') = -4j(j + l)Xj(0), 



(26) 



to write 

(V'ooh)(0) 



^X,(0)^o(j) [C\j)-C\j,)] = 



-4V2-jo(jo + l) 



M<t>) = 



(27) 



E 



Jo e 



2N sin ( 



■ sin ( d 



(0-770)) Sin 2 (0-770) 



cot 2(0 — rjO) — — sin 2(0 — rjO) + dj^ cot ( djg ^ 



For later use, let us rewrite this expression as 

-2^sin^(0-r,e) 



(^o°ii)(0)-EE 



a 



2N sin ( 



d.o(<A-';«)sin2( 



776I) [ cot 2(0 - 776I) - — sin 2(0 - t?^) + d. 
L a 



(28) 



Before proceeding, let us add a few remarks on the new boundary state, which make it particularly appealing. 

• The spin and the angle entering the state, which represent a discretised version of the conjugate intrinsic and 
extrinsic curvature variables, are conjugate variables in a precise mathematical sense: they are conjugate with 
respect to the SU(2) harmonic analysis. 

• The state is a real quantity. This is due to the fact that the phase is given by a cosine, and not a single exponent, 
as in P?)) . The fact that using the cosine and not just a single exponent does not spoil the leading order of (fTT|) 
was proved in Furthermore, one could also consider a sine term, or equivalently an SU(2) character. In 3d, 
this can be related to particle insertions [1, [lB| ■ 

• The boundary state has an interacting structure between the variables and the background; it reduces to a 
Gaussian only in the limit Jq i— > 00. For this reason it will also contribute to the corrections to the free 
propagator. In particular, the 3d analysis of this state shows that its contribution interestingly reduces the 
magnitude of the next to leading order. 
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The boundary state (j20p is factorised into contributions from single links of the boundary spin network. This 
makes its analysis much simpler than (jl7[) . but also means that it is straightforward to extend it to arbitrary 
triangulations. 



• In the definition (|20|) we took the same a for all links. Nothing prevents us from considering arbitrary configu- 
rations with a different parameter for each link. Because this parameter is in a sense related to the knowledge 
of the boundary geometry, taking the same value for all links amounts to a homogeneous description of the 
boundary state. 



V. THE 2-POINT FUNCTION AS AN INTEGRAL OVER SU(2) 



Thanks to the factorisation of the new boundary state, we can rewrite p9|) as 

1 1 f -r-r -TT 

Wab = -ijj J Y[dg,, Y[li{gi), 

where gi = g~(i)gt(i), and 



(29) 



Xji ig) V'o [jl] = ^0 (g), if l^a, b, 
E,, Xj, {gi)i'o[ji] Hii) = (V-o ° ]h)(5) if l = a, b, 



are given respectively by (|23|) and (P7)) . This expression gives the graviton propagator as an integral over SU(2). 

The choice of the adge amplitude Ae{jf) — (H/'^i/)"^ made above simplifies the computation of the two Fourier 
transforms. More generically, if Ae{jf) does not couple the spins and is simply a power of H/ dj^ , Wab still factorises 
as above, with the extra powers of dj, acting as differential operators dg in computing the Fourier transform. On the 
other hand, if the edge amplitude couples the spins and introduces some dj^j^j^j^ factors (see footnote[2|), the situation 
is different. For a single 4-simplex, these factors can always be compensated by introducing suitable counter-factors 
in the boundary state, and thus the form (j29p can be restored. For multi-simplices configurations, this will not 
work. However, for arbitrary triangulations it is not straightforward to obtain the integral representation (j29p of 
the graviton, even with the simplest choice A^^jf) = (Il/'^j/)"^- fact, the one-to-one correspondence between 
SU(2) characters and factors pT|) of the boundary state is in general lost. This makes it harder to perform the sums 
explicitly. We postpone the study of arbitrary triangulations to future work, nonetheless let us notice here that the 
correspondence is preserved by any n-valent vertex, as long as the vertex amplitude is provided by the relativistic 
spin network Therefore, an integral expression like (|29|) can be obtained for any triangulation that can be coarse 
grained to a single n-valent vertex. 

To explicitly evaluate the perturbative expansion in jo, it is convenient to use the measure ([8]) in terms of the ten 
angles. Notice that the factors sint/); in ([5]) simplify with the ones coming from '(/'o(5) and (tAo ° li)(5)j and we can 
write dMl) as 

^-fc = E E / n^<^' sUciGij) J]e,e-"^^o(^<-me) ^^^^^ ^^^^^ ^_22o (30) 

where we have absorbed the constant (2iV)^'^ in the normalisation, and we have introduced the notation 



Jo 



sin2(0a - VaS) 



cot2( 



Vat 



Jo 



sin2(0Q - ilaO) + dj„ iea 



(31) 



The expression (|30p is the contribution to the graviton propagator coming from a single 4-simplex. 

Let us stress that ([50]) is an integral over SU(2) with no sums involved, as opposed to This result is particularly 
important from a numerical point of view. The integral is in fact much easier to handle numerically than the sums, 
and the formulation provided here gives substantial progress in the numerical simulations [l^ . As far as the sum over 
the 77 signs are concerned, we show below that each configuration gives the same contribution, so that the sum gives 
a trivial redundancy that can be reabsorbed in the normalisation. 
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To define the perturbative expansion in £p wc can use tlie parameter jo ^-s above. Notice that jo enters the 
exponential of the integrand. Therefore the leading order in the jo s- oo limit can be computed using the saddle 
point approximation of the integral. We do so in the next Section. 



VI. THE SADDLE POINT APPROXIMATION 



The study of the leading order and corrections to the 2-point function formulated as a group integral requires the 
analysis the saddle point (or stationary phase) approximation of the following type of integral: 



f dii[(t)i] F{(t)i) e'^,'^irf.o(^'-me) g- Ei^^in' (01 ^ 



where e/ = ± are signs, F{(j)i) an arbitrary function (with no dependence on jo) and the measure d^,{(j)i) is defined in 
([8]) . For large jo that integral is dominated at the fixed points of the action 



2 

- s\v?{(t)i - 7]ie) - 2iei{(l)i - i]ie) 
a 



K 



i—det G/j, (32) 
jo 



where we used a Lagrange multiplier k to introduce the constraint det Gij ~ 0. To compute the equations of motion, 
notice that ^^'j}^^''' ~ —2 sin0; A*^'', where A^'^ is the determinant of the off-diagonal Z-th minor of G/j, obtained 
removing the line and column corresponding to one of the two cosi/); appearing in it. 
Then the equations of motion read: 



■ sin 2(0; -77(6') - 2ie, -2i-^ sin^jA^') VZ, 



(33) 



det G/j = 0. 
The only real solution to this set of equations is: 



0; = + 2^ '^' = (34) 

where e = ± is a global sign (independent of the link Z), and it requires the symmetry property ui = aij = criaj, with 
(7/ and aj independent signs. We have sin^; = rjiai sin 6* and A(')(0/) = ai Aq, Aq := cos6' (1 — cos 6')'^. Consequently, 
the Lagrange multiplier takes the value 

^ - ' (35) 



sin Aq 



The action evaluated on the fixed point p4[l is purely real and vanishes, S[(j>i,'K] ~ 0. If e; ^ ej, there is no real 
solution. The saddle point approximation selects two specific terms out of the sum over all possible sign assignments. 

More generally, the first part of the action, S'o[0, k] = — 2ie;0; — i-^ det G/j, is the standard action governing 
the asymptotics of the (equilateral) {lOj} symbol. It has two non-degenerate fixed points c^i = Q and 0; = —0, but 
also degenerate fixed points at (/); = 0, tt. These degenerate fixed points dominate the asymptotic behavior of the 
{lOj} symbol [U, [H, [l^ , and are responsible for the D term in (fT4|) . However, we see that their contribution to the 
2-point function is suppressed: in ()30p . there is an additional Gaussian weight, which is maximized by (/); = but not 
by = 0,7r. Therefore it kills the degenerate configurations and peaks the asymptotics around the non-degenerate 
semiclassical configuration. This provides the mechanism to neglect the D term, as done in [l|. 

We have not shown that there is no complex solution to these equations although we believe there is none.^ In 
particular, there could be fixed points close or along the imaginary line 0; — r]iQ S iM, for which the Gaussian would 



* At least, it is fairly straightforward to check that there is no purely imaginary solutions </); S ilR. 



11 



blow up. Such points can be easily found if dot G/j 7^ is allowed, but otherwise we have not been able to find any. 
The results of our partial analysis suggest that the constraint dot G/j — protects the integral from the presence of 
exponentially enhanced fixed points, a mechanism that could be relevant also for the asymptotics of the {lOj} symbol 
alone. 

Before proceeding, let us add an important remark.^ After the identification of the fixed points, one usually 
computes the integral along the complex contour such that the phase (namely the imaginary part of S) varies the 
least. This steepest descent method is implemented imposing that the imaginary part of the second derivative of S 
vanishes. However, in the following we will use the real contour, which does not satisfy this requirement (as can be 
immediately seen from the explicit expression of the second derivatives in ([55]) ). Nevertheless, as we mention at the 
end of this section, a numerical check has been performed independently from our analysis and the fit is (at least) 
of two decimals. This means that most likely there are no other (relevant) fixed points and that the real contour is 
enough for the study of the correlations. 

The fact that complex fixed points, if they do exist, are not relevant for our analysis may be due to the following 
reasons. First, notice that we are summing over all signs rji, in such a way that only the real part of the integral 
matters (see below), a fact that might lessen the relevance of the phase. Then, recall that the precise normalisation 
of the integral is not relevant but only its ratio with Af, the integral without the h insertions. Finally, we are merely 
computing the first order of the large jo asymptotics of Wab- The existence of complex fixed points might become 
relevant when computing the higher order corrections. We leave this question open for future investigations. 

Based on this preliminary analysis of the stationary points, we are now ready to state the main result of this paper 
on the asymptotics of the 4-simplex 2-point function. 

Theorem 1. The leading order of the large jo expansion of i3U\) is 

Wabijo) = (36) 
30 

where fat has only three independent entries. 

These correspond to the three independent projections of the graviton propagator along couples of normals belonging 
to the same triangles (see discussion in Section HIT)) , and should not be confused with the physical independent com- 
ponents. The explicit structure of fab{oi)^ reported in the proof below, is rather complicate, reflecting the complicate 
geometrical structure of a 4-simplex, even in the equilateral case. 

Proof. The proof goes as follows. 

1. We rewrite ([50]) introducing a Lagrange multiplier k for the constraint det G/j, 

^-"^^TrY.Y.f^^f U'^^' in A e-^°^[^-''l+'^'^'(*'-'"^) li(0a) h(0,), (37) 
■^o-'V ^,=±,,=±^ J I \ I J 

where S[(l>i, k] is given by (j32p . The normalisation A/" is given by the same quantity above without the insertions 
]h((/)a) h.{(j)b) and without the constant factor In the calculations below, we will use Af to reabsorb a number 
of overall constants, without affecting the final result. 

The expansion parameter jo enters the exponent in front of the action 5, and the field insertions h. The saddle 
point of the action is 

(pi = (f>i, <^i=Tu K = K, (38) 
where and ej are given in (|34p , and k in (|35p . 



^ We thank our referee for pointing this out. 
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2. To compute the leading order of ([37]), we expand the action to second order around the saddle point (where 
both the zeroth and the linear orders vanish), and evaluate the rest at lowest order. In particular, at the saddle 
point we have 



I saddle 



VI, 



saddle 



n 



(39) 



Consequently, we obtain the global sign J^j rji ai . When we perform the sum over the 77/ signs in (|37|) , we obtain 
identically zero unless we have the matching rji = a ai VZ, where a is an irrelevant overall sign. In the following, 
we take ai = rji. 

The leading order of ([57)1 is thus 



17!=±£=± I 



H^a+Ha)H^b + 5c^b). (40) 



3. In principle, the value of the field insertions at the saddle point is enough for the leading order. However, we 
have ]h(0) = 0, thus we need to expand the field insertions. This can be done in two steps. First, we keep only 
the terms in in (pij). thus 



ll((^a) ^{-? sin 2{(j}a - TlaO) I - siu 2(0, - -qaO) 

a L 



(41) 



Then, we expand 
write 



We have sin I (1 — Oa)''^ + "i-^ipa I — sin 254>a — 2 (50a- Using also = e ffc, we can 



ll(0,+<50a)^(-)'2 50a 

a 



2aieria-2 S(j)a 



(42) 



However, these two terms are of different orders. In fact, notice that is a Gaussian integral with width 
proportional to l/VJo- Therefore 50^ ~ 1/jo, so the second term in (|^^ can be neglected, and finally 



h{<l>,+6^a)^^ai{^-^)^eTJaScf>a. 

a 



(43) 



4. The second derivatives of the action have the following form, 



saddle 



_ 2i 

saddle JQ 



■ iji sine Aq. 



(44) 



The explicit form of A^, is given in the Appendix, and it satisfies A^^, = A^, . 

Notice that the constant factor sin6' Aq can be neglected rescaling the definition of Sk, and then reabsorbing it 
in the normalisation Af of 



5. The leading order of ([37|l is thus the following Gaussian integral, 



Wab 



-jfT^ E E^-^f- [^^'^ I n 



" 4-4- 

16 1 



; d0a (xpb exp 



-■qiriv^-^A\i,5(j)i5(l)i, -iijiScpiSK 



Z(0) dJadJt 



Z{J) 



j=o 



(45) 



where we have introduced the generating functional 



Z{J) ^ dSn Wd5(j)i exp 

m=±f-=± I 



(46) 
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6. The generating functional is a Gaussian integral that can be evaluated straightforwardly. To maintain the 
explicit symmetry of the problem, we perform first the integral over the ten angles, obtaining 



Vi=± e=± 



dSK exp 



2jo ^ 



(47) 



Here we used the fact that det yqi rjii A^, j = det A^, . Observe that consequently the sums over the rji signs give 
a trivial redundancy 2^". The remaining integral is also straightforward, and we finally obtain 

Y,{A-')wJiJv- ^ ' 



Z{J) ^ 2^° 



(27r)5 



exp 



1 



(48) 



Because A^, = A^, , the sum over the e sign amounts to taking (twice) the real part of the summand. Reabsorbing 
the irrelevant constants in the normalisation Z{0) and defining Aui = ^^/, we can then write 



Z{J) = Re 



^det Aw Ewi^-') 



exp 



1 



(49) 



7. As shown in the Appendix, the matrix {A ^)ui satisfies = fiic^) V^, where /i(a) = 2C^^+5') ' 

thus Eu'iA-^)u' = 10/i(«). 
Using this, we can finally write 



and 



AAa2jo \^det^ 10/i(a) 
8. The theorem is proved with 

fab{a) = 



AA = Re 



^det A 10/i(a) 



(50) 
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(51) 



□ 

If seen as a matrix, fabioi) has only three independent components, as it is shown in the Appendix. Geometrically, 
these are related to the three cases (i) when a and b are the same triangle, (ii) when they share a side (they are 
"adjacent" in the 4-simplex), or (iii) when they share only one point (they are "opposite" in the 4-simplex). 

For instance, choosing the value a = 1/2, we have 

W^adj (jo) ^ — , Wopp (jo) - — . (52) 
Jo Jo 

Remarkably, these values can be confirmed by a numerical analysis of (P7|) [l3 |. 
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VII. CONCLUSIONS AND OUTLOOK 

We have introduced a new boundary state to construct the graviton propagator. This new boundary state is given 
in ((20|) . and in the large spin hmit it reduces to the same Gaussian considered in [l|. Using the new state we were 
able to write the (contribution from a single 4-simplex to the) graviton propagator as an integral over SU(2). This 
is explicitly given in pop . The result is very useful from the point of view of numerical analysis (see [T3|). which is 
significantly simplified with respect to the expression of [l[. Also, the integral expression of the graviton allows to 
recast the £p expansion into the saddle point expansion of the integral. Here we evaluated the leading order, given 
in (|36p . The reason why the degenerate configurations of the Barrett-Crane model can be neglected emerges clearly: 
they correspond to non absolute minima of the action. As such, they do not enter the perturbative expansion. 

These results provide a starting point for further developments. Let us mention a few which we regard as particularly 
interesting: 

• Compute the next to leading order correction. This is a formidable task, but the integral expression here 
obtained allows a precise setting to do it. One has to consider higher orders in the expansion of the action 
P2p around the saddle point, as well as higher orders in the expansion of the field insertions. Furthermore, we 
expect different choices of edge amplitudes to affect the next to leading order (as it happens in 3d Q), and thus 
computing it should allow to discriminate between the different versions of the Barrett-Crane model: different 
choices of coupling between the 4-simplices will lead to different higher orders. From this point of view, the fact 
that the expression can also be studied numerically provides a crucial support. 

• Compute the angle correlations. This permits to reconstruct the whole tcnsorial structure of (|10|1 . and check the 
number of physical degrees of freedom. Up to now, we have looked at correlations between the spins j associated 
to two triangles of the 4-simplex. These area correlations correspond to components of the type Waabb{x,y)- 
To access the remaining components, we need to look at observables involving four different triangles of the 
4-simplex. This corresponds to correlations between dihedral angles, i.e. between two intertwiners within the 
4-simplex. Looking at such correlations was also suggested in [l7j in order to study the possibility of long-range 
correlations in BC-like spinfoam models. 

• The boundary state has a parameter, a, in some sense related to the knowledge of the boundary geometry. Can 
we consider also the angle 9 as an external parameter, to be related to the choice of triangulation? This is an 
interesting question that can be addressed looking at what happens if we choose an angle 6^6 — arccos( — 1/4) 
in the boundary state ((20)) . The phase of the boundary state would not match anymore the oscillating term of 
the {lOj} symbol and the saddle point approximation would fail. This means that the single 4-simplex would 
not give the leading order of the correlation anymore: a different spacetime triangulation is needed to obtain the 
1/jo behaviour. For example, for an angle 9 = n9 we can imagine a configuration with n 4-simplices to dominate 
over the single 4-simplex configuration, and possibly give the 1/jo behaviour. Notice that the new triangulation 
giving the leading order would depend on the group field theory coupling constant A and the choice of coupling 
between 4-simplices (edge amplitude). More work on many 4-simplices configurations, following the original 
analysis of Q , is needed in order to understand these situations and possibly highlight the precise physical role 
of A. 

• Consider the Lorentzian case. Now the relevant group is the non-compact Lorcntz group S0(3, 1) ^ SL(2, (D). 
We expect the procedure proposed in this paper to extend directly, once the harmonic analysis for the Lorcntz 
group is used. In particular, the Gaussian exp(— sin^ (j)) for SU(2) used in the boundary state should be replaced 
by the analogous exp(— sinh^ (j)) for SL(2, (D), or maybe by the simpler Gaussian exp(— 0^). The latter case could 
further simplify the analysis of the group integrals. 
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APPENDIX: GEOMETRY OF A 4-SIMPLEX 



Where more convenient, we use the double index notation for the Unks: I = IJ, where / and J are the two nodes 
hnked by I. 

Let us consider five 4d unit vectors Nj £ S^, I = 1 . . .5, and introduce the ten angles defined by their scalar 
products, cos (f)j J = Nj ■ Nj, with the convention (f>jj ~ 0. Finally, we define the 5x5 Gram matrix, G/j = cos0/j. 
These five vectors are not linearly independent, so we can find vj € such that: 



'^viNi = 0. (53) 



1=1 

This means that the 5- vector vj is a null vector for the Gram matrix Gjj. In particular, we get a constraint on the 
angles 

VJ, ^viNi-Nj^O ^ dctGij=0. (54) 

This constraint can be interpreted geometrically as follows. The five unit vectors define a unique geometric 4- 
simplex (embedded in R^) up to a global scale (4- volume of the simplex). They are the (outward) normals to the 
five tetrahedra of the 4-simplex. The closure condition of the 4-simplex reads exactly as ([55)) with the vj being the 
(oriented) 3-volumes of the tetrahedra. Furthermore, we can differentiate the equation (j54[) and contract it with the 
null vector. This gives: 

'^vi VJ sm(j)ij d(j>ij = 0. (55) 

Next, recall the well known relation Svjvj sin (pjj = iV Ajj, where V is the 4- volume of the simplex and Ajj the 
area of the triangle shared by the two tetrahedra / and J. This relation allows to write (|55|) in a simple geometric 
form as 

J^Aijd^ij^O, (56) 

which is the Schlafli identity. 

On the other hand, directly differentiating the dot Gjj = constraint gives: 

^ sin 07,7 A'-^-'M^/j = 0, (57) 

KJ 

where A*^^"^) is the off-diagonal minor obtained by removing the Ith line and Jth column from the Gram matrix G/j. 
In particular, this means that the minors are related to the areas Ajj, up a global scale factor k: 

A(") ^ n^^. (58) 
sm (pij 

This relation was used in Section II. 

The minors A*^^''^ = A^'^ play a major role in evaluating the components of the graviton propagator. In fact, 
their derivatives enter the second derivatives of the action (|32p . Using sine/); = r/iaisinO and A^'^(0j) = ctzAo, 
Aq := cos 6* (1 — cos9)^, at the saddle point ([38]) we have 

d^S 4 _ , „.K / , .rn . . . 5A«- 



0(pi 0(pv 

where 



/ -?7;6')<5o' +2i— f cos<?!)iA('UH' +sin(/), — — ) = -qi-qi, Ali, , (59) 

?n V OQl' / saddle 



= — cos 2 

saddle a ' ' ' Jo 



'4 \ 1 )9A(') 

All. ■■=[- + 2ie cot e Sw + 2ie — — — , (60) 



16 



and 



saddle 



Notice that A^^^ = . 

Evaluated at the saddle points, the derivatives of the minors 

^ ^ A 2 cos^ 0+3 cos ^+1 -r / // 

-'^i VI Ao sine cose lil^l , 



2i— sin (j)i A^'-^ ^ 2i—r]i sin9 Aq. 

7n saddle Jq 

have the foUowing values, 



(61) 



saddle 



^ ^ A l+cos^ 





Using these and the cxphcit value 9 = arccos(— i) of the angle, ([50]) reads 

if I = V, 



if / and /' adjacient, 
if I and /' opposite. 



(62) 



6 /2x/l5 i \^ 



ie ^1 if Z and I' adjacient, 
if I and I' opposite. 



(63) 



The matrix Aui := A^, has only three independent entries, coming from the three independent components of the 
matrix Consequently, each row of the matrix and of its inverse shows the same structure, which is summarised 

in the table below. For convenience, we have introduced the quantity a' 



2\/T5 
3a 



matrix A. 



II' 



/15 



along the diagonal 
occuring six times 
occuring three times 



inverse matrix {A 

y/IS (3a'+2-i)(3a' + ll8:) 

2 (3Q'+20i)(-7i+3a')(3Q'+5i) 

3VT5 (3a' + lli) 

2 (3a'+20i)(-7i+3a')(3Q'+5i) 

■ 37T5(3a' + lli) 

* (3Q'+20i)(-7i+3Q')(3Q'+5i) 



Finally, we have 



det An 



1024 
759375 



(3a' + 5i)(3a' + 20i)^{-7i + 3a')^ 



and 



Ma) :=E(^"') 
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2(2^ + 5i) 



(64) 



(65) 



from which = 10/i(«). 
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